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Abstract 

The strength of a simple soft bond under constant loading rate is studied theoretically. We 
find a scaling regime where rebinding is negligible and the rupture force of the bond scales as 
const. + (ln(kv)) 2 / 3 , where kv is the loading rate. The scaling regime is smaller for weaker bonds 
and broader for stronger bonds. For loading rate beyond the upper limit of the scaling regime, 
bond rupture is deterministic and thermal effects are negligible. For loading rate below the lower 
limit of the scaling regime, contribution from rebinding becomes important, and there is no simple 
scaling relation between rupture force and loading rate. When we extend the theory to include the 
effect of rebinding we find good agreement between theory and simulation even below the scaling 
regime. 

PACS numbers: 82.37.-j, 87.15.By 
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Non-covalent molecular bonds in biological systems are soft. Typical binding energy 
is on the order of lOfc^T. This extremely weak binding energy enables fast response to 
external stimuli which is preferred by living systems. |l| At the same time, these weak 
bonds are sensitive to thermal fluctuations, thus the study of the mechanical response of the 
bonds has to consider thermal effects. Experimentally, the strength of a soft bond is often 
characterized by the relation between the rupture force and loading rate, i.e., the dynamic 
force spectroscopy (DFS). 2] The imposed loading rate in an experiment provides a time 
scale that can be used to probe the internal dynamics of the bond under study. Therefore 
DFS experiments of soft bonds have revealed valuable structural and dynamical information 
for biological soft materials ranging from the energy landscape of several non-covalent bonds 
to the unfolding mechanisms of several types of proteins Q|. 

The theoretical problem of finding the rupture force of a bond under constant loading rate 
corresponds to evaluating the escape time of a particle in a time-varying one-dimensional 
potential well in the presence of thermal fluctuations. This is an interesting extension of 
the celebrated Kramers escape theory BQQ. To study experimental* relevant problems, 
extensions of the theory to systems with multiple bonds in parallel p] or in series p|, and 



extensions of the theory to bonds which are described by higher dimensional potentials 
have been studied recently. In the case of multiple bonds, experiments jl(| have found good 
qualitative agreement with theories [fllz]. In the case of bonds with more than one energy 
minimum, the theory [^| provides new interpretation for experimental results. 

Despite the successes on both experimental and theoretical sides, it is surprising that 



the conventional "linear theory" 
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121 ] for the DFS of a simple one-dimensional bond, 



the corner stone of theoretical analysis of most experimental studies, is challenged recently 
by two different theories that we refer to as harmonic theory 

[6] and cubic theory |13( 

respectively. Linear theory assumes that under constant loading rate the activation barrier 
for bond dissociation diminishes linear in time. The resulting rupture force scales like ln(kv) 



where ,„ is the ,oadin g rate. Recently, Hnmmer and Szabo treated the bond as a harmonic 
potential but ignores rebinding and details of potential shape near the barrier, thus we call 
this theory harmonic theory. The harmonic theory predicts a relation between the rupture 
force and the loading rate with the form const. + (Inkv ) l l 2 . At about the same time, Dudko 

n 

et. al.[13|] assumed that bond rupture occurs when the potential energy near the barrier can 
be approximated by a cubic function, thus we call this theory cubic theory. Cubic theory 
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follows an earlier work by Garg |4| ,neglects rebinding, and shows that the rupture force scales 
like const. + (In kv) 2 ^ 3 . Ref. fj showed that the harmonic theory can be extended to describe 
the unfolding of protein titin, and Ref. ljj showed that the scaling prediction in cubic theory 
for Morse potential is better than linear theory by comparing the fitt ing of rupture force 
in simulation to (ln(kv )) 2 / 3 and ln(kv) respectively. Both Ref. [f| and K| did not provide 
a quantitative analysis for the effect of rebinding. In another seemingly unrelated field of 
the physics of the atomic scale friction, linear theory ^| and cubic theory have both 
been proposed to describe the rupture of bonds (with no rebinding) between two parallel 
surfaces under relative motion. Although data from friction experiment is not able to 
discriminate between linear and cubic theory, it is believed that for atomic friction linear 
theory works for weak bonds with binding energy slightly greater than /c#T and cubic theory 
is better for bonds with realistic binding energies which is on the order of 10k bT }l5, 17|. 

To discriminate between linear, harmonic, and cubic theory, and to provide the underlying 
physical picture, including rebinding, for the rupture of molecular bonds under constant 
loading rate, in this Letter we present an analysis of simple bond DFS. Unlike Ref. flQ, we 
do not discuss higher dimensional problems. We focus on simple bonds, we consider loading 
rate beyond the scaling regime, and we consider the effect of rebinding in our theory. Our 
result shows that an extended cubic theory which takes rebinding into account is suitable to 
describe the rupture of simple bonds even at low loading rate beyond scaling regime. 

In a typical DFS experiment, a bond attached to a surface on one end and a spring on 
the other end is pulled out by the spring with a constant velocity v until the bond breaks 
and the spring recoils to its rest position. For a sufficiently simple system, a single reaction 
coordinate x, usually the distance between the particle and the substrate, is sufficient to 
describe the dynamics of the system. Let Uo(x) be the potential of the bond and Lq be 
its equilibrium position in the absence of external forces. The potential of the spring is 
U sp ring(x,t) = |(ft + L — x) 2 for t > 0, where k is the spring constant. When the total 
potential U(x,t) = Uo(x) + U spr i ng (x,t) is bistable, the minimum at smaller x is due to 
the bond, and it is denoted as x a (t), the minimum at larger x is due to the spring, and it 
is denoted as Xf(t), the barrier is denoted as x^it). x a (t), x&(£), and xj(t) are determined 
by solving — kvt = — (dUo/dx) — k(x — L ) = F(x). As Figure 1 shows, the solutions to 
—kvt = F(x) at time t are the crossing points between the stationary curve y = F(x) and 
the moving horizontal line y = —kvt which scans downward with constant rate kv . Between 
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t\ and t 2 , —kvt = F(x) has three solutions and U(x,t) is bistable. When t < t\ the bond is 
stable, when t > t 2 the bond is unstable. We further denote the local minimum (maximum) 
of F(x) as xi (x 2 ),i.e., —kvti( 2 ) = F(x 2 (i)). 

The dynamics of the system along the reaction coordinate is described by the Langevin 
equation, 

dx dU(x, t) ,, „ 

Here 7 is the damping coefficient, ( is the thermal noise with zero mean and variance 2k bT^. 
For simplicity, we assume 7 to be independent of x. The energy unit is chosen to be ksT, 
L is unit length, and unit time is chosen such that 7 = 1. When U(x,t) is bistable, the 
system is in the bound state if x(t) < Xb(t), otherwise it is in the free state. The probability 
4>(t) that the system is in the bound state at time t satisfies the reaction equation 

^ = -k off (t)<t>(t) + k m (w - m, (2) 

the probability that bond rupture occurs at time t is 

p b (t) = <P{t)k off {t)e-^ dt ' k ^ = (l - ^^'^(tOe-^^"^)^^*")^'^ k off (t)e-ti 2k °^ dt \3) 

and the mean rupture force is 

(k{vt + L -x a {t)))= [ dtP b {t)k{vt + L - x a (t)) (4) 

It is reasonable to assume that time scales to reach local equilibrium in both bound and 
free states are much smal 

Kramers escape theory (3,0] states that the on- and off-rates satisfy 



ler than time scales of transitions between the two states. Thus 
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k{t) = -L jU»{x a , t) I U"(x b , t) I e-Vto-W*-)) , (5) 
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where U"(x,t) = d 2 U(x,t)/dx 2 . For k(t) = k ff(t), x a = x a (t); for k(t) = k on , x a = Xf(t) 

Because bond rupture should occur at t closer to t 2 than ti, from Fig. 1 it is natural to 
approximate F(x) by a quadratic function near X\, and by a straight line near Xf(t). The 
resulting approximation for F(x) is 

*■(*)«< + h ' X< Xl (6) 

—k(x — Lq), x > x 2 . 



Here c = - , and /„ = - {^£) X=X1 - k( Xl - L ). Eq. © corresponds to a 

V / X=X\ 1 

potential that is cubic in x near Xi, and is the spring potential near x = L + ft. That is, 
cubic theory |l3j , not linear or harmonic theory, follows naturally from observing the generic 
shape of F(x). A straightforward calculation leads to Kramers on- and off-rates in the cubic 
theory. 



k on (t) - ^ 



l~yjk(2c(-f -kvt)) 1/2 exp |-| 0(-/o - kvt^J ' - \{x x -Lo- vtf - U {x : 
k off(t) = —\/2c(-f - kvt) exp <^ -— ( -(-/o - 
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At sufficiently large fcf rebinding is unlikely and 
= *«,//(*) x exp f / k off (t')dt' 



V~ 2c (fo + j /32/ r , ^3/2 C - A /H"(_f n -fc 1)i )3/2 1 , - 

W " 2^ exp|-^-(-/ o -fa,0 ^-_eV 9c ( ^) ^ (8) 

and the rupture force then has the simple form 
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X=X\ 

The second term on the right hand side of Eq. Q is k{x\ — x a (t)), in most cases it is 

his 



13] . 



negligible, thus cubic theory predicts the scaling relation F max ~ const. + (ln(/ct>)) 2 / 3 . 
scaling relation has been tested by data collapsing for simulating Morse potential in Ref. 
Here we provide direct comparison between simulation and theory for a couple of potentials 
so that any discrepancy between theory and simulation can be examined critically. Figure 
2 compares the rupture force predicted by Eq. © and numerical simulations for Morse 
potential U (x) = H/([l-exp(-26(x-L )/L )] 2 -l) with b = 1.5 and a power law potential 
Uq(x) = ^V{ 2 ( x ^L ) (i ~ 2(s/l ) 2 ) w hich is taken from Ref. Q, the spring constant is chosen to 
be k = 0.8. Our result shows better agreement between Eq. (J0J) and simulations for stronger 
bonds and deviation of simulation from cubic theory at high and low loading rates. 

The deviation of simulation from Eq. © at high loading rate can be understood by 
noticing that according Eq. (jSJ), the probability that bond rupture occurs between t± and t 2 
is 

P h {t)dt « 1 - e'^, (10) 
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i.e, at high loading rate with c/4nkv < 1, bond rupture is not likely to occur between t\ 
and t?,. This is because from t to t% x a (t) moves a distance X\ — x a (t) = ^J2kv(t 2 — t)/c, 
and the displacement of a particle under free diffusion during the same time interval is 



Xd(t) = \/2(t 2 — t), thus c/Ankv ~ {xd/{xi — x a (t))) 2 compares free diffusion and the motion 
of x a {t) due to the spring. When c/Ankv > 1, i.e., x\ — x a {t) < Xd(t), free diffusion near x a (t) 
is faster than the motion of x a (t), therefore the particle senses the barrier and bond rupture 
is due to barrier crossing, and Kramers theory is appropriate to describe bond rupture. 
On the other hand, if X\ — x a (t) > Xd, free diffusion alone cannot catch up the motion of 
potential minimum, therefore the particle is located near x a (t) at all time until the barrier 
vanishes. That is, thermal effect is negligible in this case. The filled symbols in the high v 
region of Figure 2 shows where c/4irkv > 2, and that is when Eq. @ starts to deviate from 
simulation result. 

The deviation of simulation from Eq. Q at low loading rate is due to the effect of 
rebinding. When the chance of observing a rebinding event between the mean rupture time 
(t r ) predicted by Eq. (jEJ) and t 2 is not negligible, i.e., 

-t 2 

k on (t)dt 
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Ajrkv 

> Oil), (11) 

rebinding has to be included in calculation for F max . Filled symbols in the low v region 
of Fig. 3 are where k on (t)dt > 0.1. For both Morse potential and power law potential, 
rebinding is important for weaker bonds. The weaker bonds not only have smaller upper limit 
for applying Kramers rate theory but also have greater lower limit for neglecting rebinding, 
therefore the deviation of Eq. Q is more serious for weaker bonds. The rupture force F max 
in the regime where rebinding is important can be calculated numerically by extending cubic 
theory to include nonzero k on in Eq. ((7j) into Eq. © and (@J). The long dashed lines in Fig. 2 
shows prediction of this extended cubic theory in regions where rebinding is not negligible. 
Indeed extended cubic theory provides a good prediction for rupture force in regions where 
Eq. (0) deviates from simulation result. 

It is interesting to notice that at low loading rate, there is still some difference between the 
rupture force in extended cubic theory and simulation. This reveals the limit of theoretical 



description for DFS of simple bonds. Figure 3 shows power law potential with W = 10, the 
case where the difference between extended cubic theory and simulation is most significant, 
and the potential predicted by cubic theory. Since the cubic theory potential is an expansion 
of Uq(x) near x±, it does not fit Uq(x) well for x smaller than L Q . Thus extended cubic theory 
does not have high accuracy when bond rupture time t predicted by extended cubic theory 
has x a (t) < Lq, and that is when 

F max <^(a; 1 )-|(x 1 -L ) 2 , (12) 

where F max is calculated from the extended cubic theory. The cross marks in Fig. 2 are where 
Fmax = ^o(^i) — 1(^1 — -^o) 2 , indeed they mark the limiting loading rates for high precision 
prediction of extended cubic theory. Finally, we note that in analyzing experimental data, 
the fitting to extended cubic theory uses data both inside and outside the high precision 
region. Therefore the resulting information like binding energy and k on should be very close 
to true values. This explains why data fitting in Ref. [l^J is so successful for Morse potential. 
Thus, extended cubic theory, an attempt to mimic the potential landscape of a bond near 
where bond rupture occurs, not only captures the basic physical picture of bond rupture, 
but also provides good quantitative result for analyzing experimental data. 

In summary, we have studied the DFS of a simple soft bond theoretically. Our analysis 
shows that rebinding is important at low loading rate, thermal effect is negligible at high 
loading rate, in the intermediate loading rate the rupture force scales as const. + (hi(/ct> )) 2 ^ 3 . 
The numerical simulations for a couple of model potentials with typical binding energy 
show that the prediction of scaling law is better for bonds with stronger binding energy. 
Rebinding is important at low loading rate, and extended cubic theory fixes the deviation 
between Eq.(IHl) and simulation. Our study shows that extended cubic theory provides the 
basic physical picture of the rupture of simple soft bonds. It would very interesting to extend 
the present model to multiple bonds in series and multiple bonds in parallel and compare 
to the analysis done in previous studies jf^O], and that is our future direction. 
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Figure captions 



• Figure 1: The solid curve is F(x). When t% < t < there are three solutions to the 
equation — (dUo/dx) — k(x — L ) = —kvt, they correspond to the potential wells and 
the barrier between them. 

• Figure 2: (a). Rupture force for Morse potential Uq(x) = W([l — exp(—2b(x — l))] 2 — 1) 
with b = 1.5 from both numerical simulation and Eq. (|9*|). (b). Rupture force for a 
power law potential U (x) = W(^-g from both numerical simulation and Eq. Q. 
Simulation result for W = 45 are presented as triangles, circles W = 30, diamonds 
W = 10. Filled symbols in the high v region are where Kramers rate theory cannot 
describe bond rupture, filled symbols in the low v region are are where rebinding is 
important. Thick dashed lines are from extended cubic theory, other lines are from 
Eq. 0. Cross marks are where F max = Uq(xi) — — L ) 2 , i.e., extended cubic 
theory does not predict F max with high precision when v is smallar than the cross 
marks. There is no cross mark for Morse potential with W = 45, in this case cubic 
theory agrees with simulation with high precision within the range of our simulation. 

• Figure 3: Solid line: power law potential with W = 10, dashed line: potential predicted 
by cubic theory. The prediction of extended cubic theory is not highly accurate when 
x a ((t r )) is located at region where cubic theory does not fit Uq{x) well, i.e., when 
x a ({t r )) < L = 1. 
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